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Abstract. In dense stellar systems, dynamical interactions between objects lead to fre- 
quent formation of exotic stellar objects, unusual binaries, and systems of higher multiplic- 
ity. They are especially important for the formation of low mass X-ray binaries (LMXBs), 
which are not only formed 100 times more efficiently than in the field, but also have a 
puzzling dependence on metallicity. In this contribution we review how compact objects 
are formed and retained, the mechanisms of dynamical formation and the specifics of the 
evolution of mass-transferring binaries with neutron stars and black holes in globular clus- 
ters - those two kinds of compact objects have different favored paths to become luminous 
in X-ray. We describe how stellar evolution affects ostensibly purely dynamical formation, 
producing the observed metallicity dependence for LMXBs. We also discuss the next puzzle 
to be solved on our journey to understand the link between LMXBs and millisecond pulsars 
formation. 
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2. Retention of compact objects 

2.1. Neutron Stars 

NSs are usually formed via collapse of the core 
of a massive star during a supernova (we will 
refer to this channel of NS formation as to a 
core-collapse, CC). The mass range in single 
stars to form a CC NSs is ~ 8 - 21M for 
solar metallicity Z=0.02 and ~ 7 - 19M for 
Z=0.001. CC channel is an NS formation path 
which is accompanied by natal kicks that have 
a very high mean velocity, ~ 400 km/s (Hobbs 
et al. 2005). Discovering natal kicks immedi- 
ately posed the problem: how NSs can be re- 
tained in GCs, where the escape velocities v esc 
are several dozen of km/s. Indeed, a "typical" 
GC of 2 x 10 5 MO, with v esc = 40km/s, will 
be able to retain only one NS if its stellar pop- 
ulation was entirely consisting of single stars, 
and only 15 if all stars were initially in binaries 
(Ivanova et al. 2008). This is in direct contra- 
diction with observations, where we know GCs 
that contain several dozen of identified radio 
MSPs (Freire 2012), and some are estimated to 
have more than a 100 (e.g., Terzan 5, Bagchi et 
al. 2011). 

An alternative formation channel of NSs is 
via electron captures which start when a de- 
generate ONeMg core reaches a critical den- 
sity - this is known as electron capture super- 
nova (ECS, Miyaji et al. 1980). ECS can occur 
during normal stellar evolution when a degen- 
erate ONeMg core is ~ 1.36M© (evolutionary- 
induced collapse, EIC). The critical condition 
can be reached in a very narrow range of initial 
masses of single stars, where the exact values 
of the mass range depend on metallicity as well 
as on underlying assumptions how mixing and 
convective overshooting should be treated, and 
hence vary between the stellar codes. Being in 
a binary helps. First of all, the mass range of 
stars where EIC can proceed is larger, due to 
mass exchanges (Podsiadlowski et al. 2004). 
Secondly, a degenerate ONeMg core with the 
critical central density could also be formed as 
a result of an accretion on a WD (accretion- 
induced collapse, AIC), or be formed as a re- 
sult of either a merger (in GCs, also due to 
a collision) of two WDs (merger-induced col- 
lapse, MIC). 



In a field population of single stars, from 
10% to 15% of all NSs ever formed can be 
formed via EIC (Ivanova et al. 2008). Fraction 
of AIC NSs and MIC NSs is rather insignif- 
icant, in case of 100% primordial binaries 
they provide ten times fewer NSs than EIC. 
However, since ECS NSs are born with much 
lower natal kicks (Kitaura et al. 2006), a sig- 
nificant fraction of them can be retained in 
GCs. The mentioned above "typical" GC can 
retain as many as 200-300 NSs, 47 Tuc cluster 
could retain up to a 1000 of NSs (Ivanova et 
al. 2008). As a result, in GCs, ECS NSs are ex- 
pected to exceed CC NSs by 30-200 times, pro- 
viding enough NSs to explain the observations. 
Since ECS NSs are expected to be at a low- 
mass range, 1.22 - 1.27M , future mesaure- 
ments of masses of NSs in GCs can verify this 
hypothesis on the prevailance of ECS forma- 
tion path for NSs in GCs. 

2.2. Black Holes 

BHs retention in GCs is less affected by na- 
tal kicks than the retention of NSs as many 
BHs are formed with kicks lowered through 
fallback (Belczynski et al. 2006). It can be es- 
timated, that, depending on metallicity, each 
15O-2OOM of stellar mass in a current "aged" 
globular cluster had produced a BH in the past 
(assuming no drastic mass loss event like tidal 
stripping occurred to that GCs). 40-50% of the 
formed BHs would be retained after their for- 
mation in a relatively massive GCS with v esc = 
50 km/s (Belczynski et al. 2006). 

It was widely anticipated in the past that 
BHs could be effectively removed by another 
mechanism - Spitzer instability would lead to a 
formation of a sub-cluster of BHs, decoupling 
BHs from the rest of GC stellar population. 
This sub-cluster is much denser than the rest of 
the GC core. Via strong interactions between 
each other, BHs quickly evaporate from the 
GCs, presumably leaving at the current epoch 
only one single or a binary pair BH (Kalogera 
et al. 2004). Detailed numerical calculations of 
such BH sub-clusters however showed that in 
massive clusters up to ~ 20% of the BHs may 
remain; and these BH sub-clusters do not reach 
equipartition (O'Leary et al. 2006). Further 
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Monte Carlo studies of a whole GC showed 
that up to 25% of initially retained BHs not 
only remained, but also participated in inter- 
actions with other stars (Downing et al. 2010). 
Most recent studies claimed that the remaining 
fraction could be up to a half of initially re- 
tained BHs which form no BH sub-cluster at 
all (Morscher et al. 2012). 

3. Dynamical formation channels 

In GCs, most of NSs or BHs lose their primor- 
dial binary companion either before, or imme- 
diately on their formation, unless it was a NS 
formed via AIC. Compact objects then can ac- 
quire a new companion via exchange into a bi- 
nary (BE, binary exchange), a physical colli- 
sion (PC) with a giant or via a tidal capture 
(TC). The formation of a binary with a NS or a 
BH does not imply that this binary will become 
an LMXB. If a formed binary can become an 
LMXB within the Hubble time while evolved 
in isolation, we call it as a direct LMXB for- 
mation. 

3.1. LMXBs with NS accretors 

TCs were found to be the least efficient mecha- 
nism of LMXBs formation in GCs; when suc- 
ceed, they form predominantly LMXBs with 
MS donors (this and further statistics in this 
subsection is taken from simulations described 
in Ivanova et al. 2008). PCs lead to a forma- 
tion of NS-WD binaries (Ivanova et al. 2005), 
which, if become LMXBs, are known as ultra- 
compact X-ray binaries (UCXBs). Most of 
TCs and PCs result in a direct NS-LMXB for- 
mation. NSs which were formed via AIC, in 
most cases, would become LMXBs soon af- 
ter NS formation as their orbit only slightly 
expanded; donors in such LMXBs remain the 
same as at the moment of AIC and usually are 
MSs. 

BEs can form a binary with a NS and any 
type of a companion, however they usually do 
not lead to a direct formation of NS-LMXBs. 
To become an LMXB, a post-exchange binary 
would have to be altered further, via subse- 
quent dynamical encounters. This usually oc- 
curs as a result of the process we term as "ec- 



centricity pumping". Eccentricity pumping is 
a result of a sequence of fly-by encounters 
which gradually boost eccentricity. A fly-by 
encounter is an encounter that occur at a dis- 
tance large enough so it does not affect signif- 
icantly a binary separation during one individ- 
ual encounter, neither is leads to a companion 
exchange. The fly-by encounters are generally 
known to gradually harden a binary, decreas- 
ing its orbital separation, if that binary was ini- 
tially hard (Heggie 1975). However, for bina- 
ries with NS, the more important role in a suc- 
cessful NS-LMXB formation was identified to 
be due to the eccentricity pumping. 

BEs predominantly lead to a formation of 
LMXBs with MS and RG donors. Formation 
of UCXBs via BEs, while less productive than 
LMXBs formation with MS and RG donors, 
is still non-negligible, and produces just few 
times less UCXBs that PCs. We note that post- 
exchange NS binaries, as well as post-PC NS 
binaries, are eccentric at the start of RLOF. 

NS binaries can also form hierarchically 
stable triples through encounters with binaries. 
E.g., in a cluster similar to 47 Tuc, each NS 
binary has a 5% chance to form a triple in a 
Gyr (Ivanova 2008). A third of these triples is 
affected by a Kozai mechanism. Kozai mecha- 
nism, if coupled with tidal friction, drives the 
inner binary to merge or RLOF before next 
interaction (Eggleton & Kisseleva-Eggleton 
2006). In numerical simulations, about a half 
of binary MSPs (bMSPs) were in triples at 
some point in their past! 

3.2. LMXBs with BH accretors 

Observationally implied formation rate of BH- 
LMXBs is very low: indeed, only one po- 
tential BH-LMXB is present at the time per 
1 - 3 x 1O 9 M in GCs (Ivanova et al. 2010). 
Hence theoretical studies of the formation of 
BH-LMXBs using a self-consistent dynami- 
cal code to model a complete GC are compu- 
tationally unfeasible. To obtain statistics and 
rates, Ivanova et al. (2010) used a simplified 
approach, where the formation of BH-LMXBs 
was studied as per a BH via separate forma- 
tion channels. The most reliably identified at 
the time GC BH-LMXB, had low Ha/[OIII] 
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ratio in its spectra and strong, broad OIII lines 
(Zepf et al. 2008), signaling about likely WD 
donor. Hence only the formation of BH-WD 
LMXBs was studied in details. 

It is hard to have a direct formation of BH- 
LMXBs: most of PCs lead to a formation of too 
wide binaries(Ivanova et al. 2010), AIC forma- 
tion is not applicable, BEs form too wide sys- 
tem (Kalogera et al. 2004). Indeed, Ivanova et 
al. (2010) found that the formation of a BH- 
LMXB is necessarily a multi-step one. First 
step is the formation of a "seed" BH-WD bi- 
nary. This occurs either via BE or a PC. Such a 
binary, with the rare exception of a very small 
fraction post-PC BH-WD binaries, can not be- 
come an LMXBs if evolved in isolation. 

Post-PC BH-WD binaries however can be 
transformed into future BH-LMXBs via triple- 
induced mass transfer (TIMT). For BH-WD bi- 
naries, the rate of formation of an hierarchi- 
cally stable triple in encounters with binaries 
was found to exceed the rate of strong encoun- 
ters with single stars, if the binary fraction is 
about a few per cent. A fraction of these hier- 
archically stable triples would be formed with 
high enough inclination to be affected by a 
Kozai mechanism. For the Kozai BH triples, 
a time-scale for reaching the maximum eccen- 
tricity was found to be shorter than a time- 
scale for the next strong dynamical encounter. 
Hence, once formed, a seed BH-WD binary 
with a binary separation < 80/? G can be quickly 
brought to the regime where gravitational wave 
radiation can shrink the binary to RLOF within 
the Hubble time. 

Most of post-BE BH-WD binaries are how- 
ever too wide to be trasformated by TIMT. 
They need to be hardened, and this usually 
occurs as a result of many encounters; evo- 
lution though a hardening sequence can take 
up to a Gyr. Most of initially formed post-BE 
BH-WD binaries get destroyed along this se- 
quence, but a few per cent of them survive 
until their binary separation is short enough 
for TIMT. (Hardening) -> TIMT gravitation 
wave radiation — > UCXB sequence was shown 
to provide enough BH LMXBs to explain the 
observed number of BH-LMXBs (this required 
that at least 1% of initially formed BHs is re- 
tained in a GC). It is remarkable that recent ob- 



servations showed that one of GC BH-LMXB 
is likely in a triple system (Maccarone et al. 
2010). 

Recently, two accreting BH-candidate bi- 
naries were detected in radio in M22 (Strader 
et al. 2012). The potential companions for both 
of the objects, based on the identified optical 
counterparts, are either low-mass MS stars, or 
WDs with the luminosities below the detection 
limit of the used HST data. The formation of 
BH-MS LMXBs hence might be re-accessed in 
a future, with more thorough numerical studies 
that employ a self-consistent method to emu- 
late a whole GCs, and possibly with a better 
understanding of TCs by BHs. 

4. Metallicity dependence 

A metal-rich GC is ~ 3 times more likely 
to contain an LMXB than a metal-poor GC 
(Grindlay 1993). This holds for Milky Way 
GCs as well as for extragalactic GCs. For ex- 
tragalactic sample this ratio also was found to 
hold across the range of X-ray luminosities, 
from 2 x 10 37 erg s -1 to 5 x 10 38 erg s -1 (Kim 
et al. 2012). For even higher X-ray luminosity, 
where the accretors are likely BHs, this ratio is 
less certain, but still is above one, 2.5^'j. 

Relatively low X-ray luminosities, L x ^ 
2 x 10 37 erg s~\ can be provided by LMXBs 
with a variety of donors, however LMXBs with 
L x <: 2 x 10 37 erg s _1 are unlikely to be MS- 
LMXBs (Fragos et al. 2009). Indeed, in Milky 
Way GCs, only UCXBs are present at L x S 
2 x 10 37 erg s _1 . Hence despite having the same 
ratio for Milky Way GCs and for extragalactic 
samples, the observed metallicity dependence 
must be analyzed separately - different L x im- 
plies different donors! 

Furthermore, in Milky Way GCs, there is 
no identified bright (L x > 10 36 erg s 1 ) LMXB 
with a MS donor in metal-poor GCs, while 
there are four bright MS-LMXBs in metal-rich 
GCs. This led to the idea that the formation 
of MS-LMXBs in metal-poor GCs is halted 
Ivanova (2006). In the MS-LMXBs the most 
effective mechanism of angular momentum re- 
moval is the magnetic braking, which operates 
only in stars with an outer convective zone. 
Metal-poor MS stars lack an outer convective 
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zone in a large mass range of low-mass MS 
stars. If formed, a metal-poor MS -LMXBs is 
expected to be in quiescence for the most du- 
ration of the mass transfer. 

The mechanism that prevents the formation 
of NS-MS LMXBs, however, can not explain 
the observed ratio of brighter, Lx ^ 2 x 10 37 
erg s _1 , LMXBs in extragalactic GCs - there 
the potential donors are likely RGs and WDs. 
Formation of NS-RG donors is predominantly 
a direct LMXB formation via BEs, where the 
evolutionary expansion of a giant brings the 
system to the start of RLOF. The formation 
rate of such systems depends on dynamical 
properties of a GC (also known as the stellar 
encounter rate), as well as on the population 
of binaries with RGs and the number of NSs. 
Formation of NS-UCXBs in GCs is predomi- 
nantly a direct LMXB formation via PC with a 
RG; post-PC binary shrinkage is solely due to 
angular momentum loss via gravitational wave 
radiation. The formation rate of UCXBs de- 
pends as well only a GC dynamical proper- 
ties and the number of encounter participants, 
which are NSs and RGs. The stellar encounter 
rate was found to be statistically identical in 
metal-rich and metal-poor GC samples (Kim 
et al. 2012). From our current understanding 
of NS formation and retention, it is expected 
that slightly more of NSs can be contained in a 
metal-poor GC than in a similar metal-rich GC 
- the opposite trend to the metallicity depen- 
dence. Hence, the only remaining quantity that 
might be different is related to RGs. 

Ivanova et al. (2012) compared how the 
populations of RGs can differ between metal- 
rich and metal-poor GCs. It has been always 
appreciated that GCs of different metallicities 
would have different turn-off masses for their 
main sequences: metal-poor MS stars evolve 
faster and hence metal-poor GCs have stars 
of lower mass at their MS turn-off. It ap- 
peared, that metal-poor stars live through RGs 
stage faster as well. Since RGs evolve faster, a 
metal-poor GCs would contain RGs that have 
a smaller mass range than RGs in a metal- 
rich GCs. Detailed evolutionary calculations 
showed that the number of stars that a metal- 
poor GCs can have in a RGs stage, assuming 
that GCs of different metallicities had the same 



initial mass function, is about 2 times smaller! 
The additional effects are: i) the lifetime of 
a RG-NS mass transferring system, which is 
shorter in metal-poor GCs; ii) the appearance 
of a RG-NS LMXBS - metal-poor systems 
are predominantly transient while metal-rich 
are more likely to be persistent and iii) the 
mass of a RG, which is larger in a metal-poor 
RGs, boosts the collision rate with RGs and de- 
creases the post-PC separations, resulting a di- 
rect LMXB formation for a larger number of 
encounters. Combined, all these effects result 
in the about 3 times more frequent formation 
of UCXBs in metal-rich clusters than in metal- 
poor. 

5. Production of radio MSPs 

The anticipated ultimate fate of an NS-LMXB 
is to exhaust available donor's material, spin- 
up an accretor and, depending on how much 
mass was accreted, possibly appear as a MSP. 
The comparison of the population of binary 
MSPs (orbital periods, companion masses, ec- 
centricities) is therefor a crucial test for LMXB 
evolutionary theory. As GCs are known to be 
LMXB factories, it has been expected that they 
would also contain a large number of MSPs. 
Extensive dedicated searches were performed 
specifically for GCs, and as a result most of 
identified anywhere radio MSPs (144 to-the- 
date) are located in GCs (Freire 2012), and 2/3 
of them are in binaries. 

The comparison of bMSPs populations 
in simulations of GCs and in observations 
showed nice agreement for the most of the for- 
mation channels: for BE (with a MSP being 
exchanged into a binary) as well as for post- 
mass transfer system with initial RG and MS 
donors. However, a striking disagreement was 
found for the theoretically expected progenies 
of UCXBs. Both simplified estimates and de- 
tailed numerical simulations of GCs were able 
to re-produce a number of UCXBs compara- 
ble to the observed number of UCXBs in GCs 
of different dynamical properties, however this 
did not occur for post-UCXBs bMSPs. 

The lifetime of an UCXB at L x > 10 36 
erg s _1 is only ~ 10 8 years. If the forma- 
tion rate of UCXBs in the past was not drasti- 
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cally different from the current, then each cur- 
rently present UCXB with L x > 10 36 erg s" 1 
implies a hundred of UCXBS formed in the 
past; those UCXBs should have become bM- 
SPs by now. Milky Way GCs contain about a 
dozen of such bright UCXBs at present, im- 
plying that there should be a thousand of post- 
UCXBs bMSPs. However, they are completely 
lacking in the observed sample of radio bM- 
SPs in GCs, both by the number and by or- 
bital period-donor mass distribution (Ivanova 
et al. 2008). A similar lack of post-UCXBs ra- 
dio bMSPs was identified at the same time for a 
field population of bMSPs (Deloye 2008). We 
note that at least one post-UCXBs bMSP was 
detected in X-rays and y-rays. (Krimm et al. 
2007). The reason for post-UCXBs MSPs to be 
radio-quiet is not yet understood and is a new 
puzzle posed by GCs for us. 

6. Conclusions 

• Population of NSs in GCs relies upon elec- 
tron capture supernovae: retained ECS NSs 
are expected to exceed CC NSs by far. 

• Detections of BH-LMXB in GCs indi- 
cated that BH population in GCs is non- 
negligible. This agrees with the results of 
most recent dynamical simulations of GCs. 

• LMXBs in different ranges of X-ray lumi- 
nosity are formed via different dynamical 
formation channels. The metallicity depen- 
dance is mainly determined by red giants 
life-times for Lx <: 2 x 10 3 7 erg s _1 , but is 
also influenced by magnetic braking for a 
lower-range of bright LMXBs. 

• Dynamical formation of hierarchically sta- 
ble triples is a crucial process for LMXBs 
in GCs with both BHs and NSs accretors. 

• Mass transfer in UCXBs does not result in 
the formation of radio bMSPs. 
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